Published online 29 April 2014 



Nucleic Acids Research, 2014, Vol 42, No. 10 6753-6761 

doi: 10.1093lnarlgku262 



A family of synthetic riboswitches adopts a kinetic 
trapping mechanism 

Dennis M. Mishler and Justin P. Gallivan' 

Department of Chemistry, Emory University, 1515 Dickey Drive, Atlanta, GA 30322, USA 

Received October 13, 2013; Revised IVIarch 17, 2014; Accepted IVIarch 18, 2014 



ABSTRACT 

Riboswitches are sequences of RNA that control 
gene expression via RNA-ligand interactions, with- 
out the need for accessory proteins. Riboswitches 
consist of an aptamer that recognizes the ligand and 
an expression platform that couples ligand binding 
to a change in gene expression. Using in vitro selec- 
tion, it is possible to screen large (-^1 0^^ members) li- 
braries of RNA sequences to discover new aptamers. 
However, limitations in bacterial transformation ef- 
ficiency make screening such large libraries for ri- 
boswitch function in intact cells impractical. Here 
we show that synthetic riboswitches function in an 
E. coii S30 extract in a manner similar to how they 
function in intact E. coii cells. We discovered that, 
although this family of riboswitches regulates the 
initiation of protein translation, the fate of whether 
an RNA message is translated is determined during 
transcription. Thus, ligand binding does not bias a 
population of rapidly equilibrating RNA structures, 
but rather, co-transcriptional ligand binding kineti- 
cally traps the RNA in a conformation that supports 
efficient translation. In addition to providing new in- 
sights into the mechanisms of action of a family 
of synthetic riboswitches, our experiments suggest 
that it may be possible to perform selections for novel 
synthetic riboswitches in an in vitro system. 

INTRODUCTION 

The ability to precisely control bacterial gene expression is 
critical for experiments in both basic and applied microbi- 
ology. Synthetic riboswitches are attractive tools for con- 
trolhng gene expression because the level of expression is 
regulated in a ligand-dependent fashion (1,2). Riboswitches 
consist of an aptamer domain that recognizes the ligand, 
and an expression platform that modulates gene expres- 
sion. Naturally occurring riboswitches are common in bac- 
teria and regulate the synthesis of various metabohtes (3,4). 



The ligands recognized by natural aptamer domains include 
amino acids, nucleotide derivatives, coenzymes, and even 
ions (3). In vitro selection techniques make it possible to 
screen large libraries of RNA (>10^^ members) to discover 
aptamers that bind to novel ligands both tightly and selec- 
tively (5,6). We and others have shown that it is possible to 
convert known aptamers, such as the theophylline aptamer, 
into functional riboswitches through high-throughput in 
vivo screening or rational design (7-13). However, both in 
vivo screening and rational design have limitations. 

Rational design of riboswitches is often based on meth- 
ods for computing the relative thermodynamic stabilities of 
competing RNA structures (14-17). While a variety of com- 
putational methods, such as mfold (18) and Vienna RNA 
(19), can reasonably predict the relative ground-state free 
energies of different RNA conformations, they are less able 
to predict the structures and free energies of the intermedi- 
ates and transition states that separate the conformations. If 
the kinetic barriers between conformations are sufficiently 
low, thermodynamic (equilibrium) models for riboswitch 
function should work well. However, if there are significant 
kinetic barriers between riboswitch conformations, equilib- 
rium models that fail to account for these barriers will hkely 
fail. To determine if these computational predictions are ac- 
curate, the predicted riboswitches must be tested in a biolog- 
ically relevant system, a few sequences at a time. 

Synthetic riboswitches have also been discovered using 
high-throughput screens or selections in cells (11,20-22). 
Such methods are attractive because they neither require 
prior structural knowledge of the aptamer, nor presup- 
pose a mechanism of action for the resulting riboswitch(es). 
While the best genetic selections performed in Escherichia 
coii can assay library sizes on the order of 10^ members 
(21,23), a limitation imposed by the transformation effi- 
ciency, such libraries are considerably smaller than those 
that may be assayed in a cell-free system. By compari- 
son, most aptamers are discovered using in vitro selections, 
which are capable of assaying libraries on the order of 10^^ 
members (5,6). 

Because cell-free systems can sample much larger library 
sizes than cell-based systems (>10^^ versus 10^ members). 



To whom correspondence should be addressed. Tel: +1 404 712 2171; Fax: +1 404 727 6586; Email: justin.gallivan@emory.edu 
Present address: 

Dennis M. Mishler, Department of Molecular Biosciences, University of Texas at Austin, Austin, Texas, USA. 
© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.Org/licenses/by-nc/3.0/), which 
permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact 
journals.permissions@oup.com 



6754 Nucleic Acids Research, 2014, Vol 42, No. 10 



it is tempting to think that it may be possible to directly se- 
lect for synthetic riboswitches in vitro (24), much in the way 
that aptamer selections are performed today. However, since 
many apphcations of riboswitches will ultimately require in- 
tact cells, it is also important to establish that riboswitches 
function similarly both in vitro and in cells (25). 

Here we show that riboswitches expressed in E. coli S30 
extract function in a manner similar to the way they func- 
tion in E. coli cells. To our surprise, it was critical that tran- 
scription was active for the vast majority of riboswitches to 
function well in extract. Our results indicate that, at least in 
extract, many of these riboswitches do not act as rapidly 
equihbrating systems, but rather as kinetic traps that re- 
quire the presence of ligand during transcription (Figure 1). 
These findings highlight the utility of cell-free extracts in the 
design and discovery of new synthetic riboswitches and the 
characterization of known riboswitches. 

MATERIALS AND METHODS 

General considerations 

All plasmid manipulations utilized standard cloning tech- 
niques and all constructs were verified by DNA sequencing. 
Purification of plasmid DNA, PGR products and enzyme 
digestions were performed using kits from Qiagen. Oligonu- 
cleotides were purchased from IDT. 

Cloning 

The luciferase constructs derived from the switch pack se- 
quences (26) were created using standard cloning proce- 
dures. The specific riboswitch sequences can be found in 
Supplementary Figures S5 and S6. The firefly luciferase 
gene was derived from the plasmid pBESTluc (Promega) 
using forward primer DMM141 (see below for sequence) 
and reverse primer DMM157. Overlap PGR was then per- 
formed using this PGR product and PGR products derived 
from the 5^ UTR sequences in Supplementary Figure S5 as 
templates. The primers for the overlap reactions were for- 
ward primer DMM161, which contains the tac promoter 
sequence, and reverse primer DMM157. 

PGR products were purified using a PGR purification kit 
(Qiagen) and then digested with BamH I and Pst I, under- 
lined below, and gel purified using a gel extraction kit (Qia- 
gen). Digested product was ligated into a pUG19 vector, fol- 
lowed by butanol precipitation and resuspension with 20 |xl 
water. TOP 10 cells were transformed with the ligation prod- 
ucts via electroporation, plated on LB agar with ampicillin, 
and grown overnight at 37°G. Golonies were picked the next 
day and grown in LB liquid media with ampicillin overnight 
at 37°G. Plasmid DNA was isolated from these overnight 
cultures via miniprep kit (Qiagen). DNA sequencing con- 
firmed the identities of these plasmids. All cloning enzymes 
used above were obtained from NEB. Luciferase reporters 
for the D family riboswitches (Supplementary Figure S6) 
were created using the same procedures and primers. 

DNA templates for in vitro transcription were generated 
from the above plasmids using forward primer DMM107, 
which contains the T7 promoter sequence, and reverse 
primer DMM157. 

Primers mentioned above are as follows: 



DMM107: 5^-TAATAGGAGTGAGTATAGGGAGTG 
AGTATAGGTAGG-3^ 

DMM 14 1 : 5^- ATGGAAGAGGGGAAAAAGATAA 
AGAAAGGGGGGG-3^ 

DMM 157: 5^-AGGATCTGCAGGGGGAGAGGAGT 
AAGGTGTGGGGTTAGAATTTGG-3^ 

DMM 161: 5^- AGTGAGGGATCCGAGGTGTTGAGA 
ATTAATGATGGGGTGGTATAATGTGTGGG-3^ . 

Luciferase assays 

In vivo luciferase reactions: saturated overnight cultures of 
E. coli TOP 10 cells containing a pUG19 plasmid with a lu- 
ciferase reporter under the control of a riboswitch were di- 
luted 1 :500 in fresh LB/ AMP media and allowed to grow to 
mid-exponential phase. Gulture was then prepared for lu- 
ciferase measurements by placing the reactions on dry ice 
and lysing with Ix GGLR (Promega) following the manu- 
facturer's procedure. 

Run-off in vitro transcription reactions were performed 
using the AmpliScribe T7-Flash kit (Epicentre) per the 
manufacturer's procedure with ~1 |xg of DNA template. 
After 30 min at 37°G, reactions were treated with DNase 
I for 15 min and then precipitated with 5 M NH4OAC, as 
described in the AmpliScribe protocol. 

In vitro transcription and/or translation reactions were 
performed using S30 E. coli extract (Promega, item num- 
bers LI 020 and LI 030). Reaction conditions followed those 
recommended by the manufacturer, and included S30 pre- 
mix, S30 extract, amino acids, template DNA or RNA, and 
either theophylline or water. An amount of 40 ng/ixl of 
plasmid DNA template or 800 ng/ixl of in vitro transcribed 
RNA template was used for each reaction. Theophylline 
concentrations were typically 2 mM, except for the dose re- 
sponse curves. Reaction volumes varied from 10 to 50 |jl1, 
but results were reproducible regardless of total volume. Re- 
actions were incubated at 37°G for 30 min, unless otherwise 
noted. Reactions were quenched by placing them on ice with 
dilution buffer provided with the S30 extract per the manu- 
facturer's instructions. 

Luciferase assays were conducted using samples from in 
vitro translation or cell cultures. In each case, ^20-30 |jl1 
of sample was then added to a well on a black opaque 96- 
well plate. An equivalent amount of luciferase assay reagent 
(Promega) was added to each well and mixed. Luciferase 
activity was detected by reading the 96-well plates with a 
Biotek Synergy HT plate reader, using the luminescence 
procedure with a sensitivity of 150 for most experiments. 

Presentation of the luminescence data: In vivo luciferase 
expression data were initially normalized relative to OD600 
values at the time of harvesting. For in vitro translation re- 
actions, the data were normalized relative to the amount of 
DNA or RNA template used. These values were then nor- 
malized relative to a single control reaction or time point 
within that experiment. More specifically, the values of each 
reaction or time point were divided by the value of the sin- 
gle control reaction or time point, generating a normal- 
ized scale of 0-100%. Normalized expression data were then 
averaged between at least three independent experiments. 
Standard deviations were determined and graphed as error 
bars. From these average expression values, estimated acti- 
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equilibrium model 




kinetic trapping model 




Figure 1. Models for riboswitch function. In the equilibrium model (top), a newly transcribed RNA can adopt a conformation in which the ribosome 
binding site (RBS) is paired and translation is inefficient. This conformation is in equilibrium with a conformation in which the RBS is free (efficient 
translation). The presence of theophylline (red star) shifts the equilibrium to favor the RBS-free conformation. In the kinetic trapping model (bottom), in 
the absence of theophylline, the nascent RNA adopts a conformation in which the RBS is paired and translation is inefficient. If theophylline is present 
during transcription, the theophylline-bound RNA adopts a conformation in which the RBS is free and translation is efficient. In the kinetic trapping 
model, the free and bound structures do not equilibrate on the timescale of the experiment. 



vation ratios (ARs) were determined and presented. Specif- 
ically, an AR is calculated by dividing the amount of gene 
expression in the presence of theophylline by the amount of 
gene expression in the absence of theophylline. 

Dose response curves were conducted in vivo or in S30 ex- 
tract as described above, but with varying amounts of theo- 
phylline as shown in the figures. 



Secondary structure models 

Supplementary Figures S5 and S6 and the estimated AG 
for riboswitch sequences were generated using mfold (18), 
specifically the Quikfold option for RNA, using the entire 
5^ UTR and the first two codons of the luciferase reporter. 
Inclusion of additional codons did not appear to disrupt the 
structures shown here. 



Decoupled transcription /translation reactions in S30 extract R ES U LTS 



In vitro transcription/translation reactions in S30 extract 
with rifampicin and 40 ng/ixl of DNA template were con- 
ducted as described above with the following exceptions. 
Amino acids were omitted during the initial 5 min in all 
cases. After 5 min (the transcription-only phase), a solution 
of rifampicin in water (2.5 mg/ml) was added, yielding a 
final reaction concentration of 250 |xg/ml rifampicin. One 
minute later, amino acids were added to the reaction as well 
as additional water and/ or theophylline to reach a final con- 
centration of either no theophylline or 2 mM theophylline. 
The reaction (now in the translation-only phase) proceeded 
for the amount of time indicated followed by quenching as 
described above for in vitro luciferase assays. 



We first tested whether a set of previously reported 
theophylline-responsive synthetic riboswitches could con- 
trol the expression of a luciferase reporter gene in intact E. 
coli cehs (26,27). These riboswitches were selected because 
they were previously published and are frequently requested 
by other laboratories (26,27). We created six plasmids (A-F) 
in which the tac promoter was used to drive the expression 
of a luciferase reporter gene under riboswitch control. E. 
coli cells were transformed with one of the plasmids, grown 
on selective solid media, and then in selective liquid media 
either in the presence or absence of theophylline. As shown 
in Figure 2A, cells harboring the riboswitch-controlled lu- 
ciferase reporter exhibit theophylhne-dependent increases 
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Figure 2. (A) Performance of riboswitches expressed in intact E. coli cells. 
Left axis: luciferase activity in the absence (open circles) or presence (filled 
circles) of 2 mM theophylline. Values are normalized to the C control in 
the absence of theophylline. The 'A control' and 'C control' have identical 
ribosome binding site sequences to riboswitches A and C, respectively, but 
lack the theophylline-binding aptamer. Error bars reflect the standard de- 
viation for three independent determinations. Right axis: the activation ra- 
tio, which is expressed as the relative luciferase activity in the presence of 2 
mM theophylline divided by the activity in the absence of theophylline. (B) 
Performance of riboswitches expressed in E. coli S30 extract (RNA tem- 
plates). Note that the right y-axis range has been changed by an order of 
magnitude to better reveal the smaller activation ratios. (C) Performance of 
riboswitches expressed in E. coli S30 extract (DNA templates). Note that 
the right y-axis range has been changed, relative to panel (A), by a factor 
of 2. Data for (B) and (C) are presented as in (A). 



in luciferase activity in the presence of 2 mM theophylline. 
Cells containing the luciferase reporter alone show mini- 
mal theophylline-dependent changes in luciferase activity. 
The behavior of these riboswitches, including their ARs, 
with the luciferase reporter gene mimics that previously ob- 
served with the p-galactosidase reporter gene (26), although 
the E switch does exhibit differences that may result from 
sequence-specific interactions between the riboswitch and 
one or both of the reporter genes. These data indicate that 
our series of riboswitch-controlled luciferase reporters reg- 
ulate gene expression in a theophylline-dependent manner 
in vivo. 

Previous results from our lab suggested that related ri- 
boswitches functioned via a post-transcriptional mecha- 
nism, in which theophylline binding biases equilibrating 
RNA structures toward conformations that favor transla- 
tion (Figure 1) (22). However, we were unable to fully rule 
out competing mechanisms of riboswitch function. Given 
these findings, we anticipated that mRNA templates pro- 
duced by in vitro transcription would yield theophylline- 
dependent changes in luciferase activity when added to an 
in vitro translation system. 

To test this, we placed the DNA constructs described 
above under control of a T7 promoter to drive the expres- 
sion of the riboswitch-luciferase reporter gene. mRNA was 
produced by in vitro transcription, and the purified tran- 
scripts were used as templates for in vitro translation reac- 
tions in E. coli S30 extract. Translation reactions were per- 
formed either in the absence or presence of 2 mM theo- 
phylline (Figure 2B). Surprisingly, all of the riboswitches 
showed little-to-no theophylline-dependent increases in lu- 
ciferase activity. Only the E switch with the luciferase re- 
porter gene showed an AR (~ 3 -fold) similar to its in vivo 
AR (~ 6-fold). In contrast, a control mRNA lacking a ri- 
boswitch showed robust luciferase activity both in the pres- 
ence and absence of theophylline, demonstrating that the 
S30 system is capable of efficiently translating mRNA pro- 
duced by in vitro transcription (Figure 2B). The quality of 
the mRNA templates was confirmed via agarose gel (Sup- 
plementary Figure SI). mRNA templates were also heated 
to 65° C in the presence of the ligand and then cooled prior 
to in vitro translation in extract. Surprisingly, we observed 
no ligand dependent changes in luciferase activity (Supple- 
mentary Figure S2). 

Since in vitro transcribed riboswitch mRNA templates 
functioned poorly in S30 extract, we asked whether DNA 
templates actively transcribed to mRNA in S30 extract 
would function as riboswitches. To test this, we introduced 
the same plasmid DNA previously used in E. coli cells into 
S30 extract that is competent for both transcription and 
translation. As depicted in Figure 2C, DNA templates pro- 
duce theophylhne-dependent increases in luciferase activity 
in S30 extract. This is in stark contrast to the results seen 
when mRNA templates were used. Although the ARs ob- 
served in extract are approximately 2-fold lower than those 
seen in cells, the riboswitches with the highest ARs retain 
this distinction both in vivo and in extract (compare Figure 
2A-C). 

Our data suggest that active transcription is impor- 
tant for riboswitch function in S30 extract, while our pre- 
vious studies suggested that active transcription is not 
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essential for riboswitch function in cells (22). To further ex- 
plore the relative contributions of transcription and trans- 
lation in extract, we decoupled the two processes. This 
was accompHshed by adding a DNA template to the 
transcription/ translation extract either in the absence of 
amino acids (Figure 3 A, first bar), or in the presence of 
a transcription inhibitor, rifampicin (Figure 3 A, last bar). 
Only when both transcription and translation can occur do 
we see luciferase activity. If either transcription or transla- 
tion is inhibited, through omission of amino acids or by ad- 
dition of rifampicin, no luciferase activity is observed. In 
Figure 3B, we present a schematic of the decoupling pro- 
cedure we use: after 5 min of reaction time in the absence 
of amino acids, during which RNA is transcribed but not 
translated, we added the transcription inhibitor rifampicin 
(22). After waiting for one additional minute, we added a 
mixture of the 20 canonical amino acids to allow transla- 
tion, while transcription initiation is blocked. Control lu- 
ciferase activity generated using this procedure is presented 
in Figure 3C. Theophylline was added either at the start of 
transcription, at the start of translation, or not at all, allow- 
ing us to determine the importance of ligand either during 
or after transcription. The results of the time course exper- 
iments are shown in Figure 3D. 

As expected, the reaction lacking theophylHne shows very 
little increase in luciferase activity over time. To our sur- 
prise, the reaction that had theophylline present only dur- 
ing translation also shows Httle increase in luciferase ac- 
tivity. This stands in contrast to the reaction in which 
theophylline was present during both transcription and 
translation. In this case, we observe robust theophylline- 
dependent increases in luciferase activity. Even after ac- 
counting for the extra time that theophylline is present in the 
transcription/ translation reaction (11 min) versus transla- 
tion only (5 min), there is a significant increase in signal that 
results from having theophylline present while the mRNA 
is being transcribed. To confirm that this difference was not 
simply a result of the time difference described above, we 
extended the translation period for a total of 30 min and 
found identical behavior (Supplementary Figure S3). These 
results suggest that to achieve the highest ligand-dependent 
changes in gene expression in this in vitro system, it is critical 
to have hgand present during transcription (Figure 3D and 
Supplementary Figure S4). However, our previous mecha- 
nistic studies carried out in intact cells (22) suggested that 
there is also a post-transcriptional component to riboswitch 
function. 

To reconcile these observations, we hypothesized that 
these synthetic riboswitches can function through two com- 
peting mechanisms: (i) a co-transcriptional kinetic-trapping 
mechanism, in which it is critical to have ligand present 
while transcription is active, and (ii) an equilibrium mech- 
anism in which mRNA structures equilibrate rapidly. We 
speculated that riboswitches with higher background lev- 
els of protein expression would have a greater contribution 
from the equilibrium mechanism. To test this, we studied 
a series of closely related riboswitches derived from the D 
switch. These sequences are identical with the exception of 
an 8 nt sequence between the aptamer and the start codon 
(Table 1). 



The predicted secondary structures and free energies of 
the OFF states (Supplementary Figure S6) were generated 
using mfold (18). We tested these switches in growing E. 
coli cells, and then plotted the observed AR against the es- 
timated — AG of the OFF state sequester helix (SH) for 
each riboswitch (Figure 4). The resulting plot shows a trend: 
as the predicted stability of the OFF state SH increases, 
the observed in vivo AR also increases. This result aligned 
with our expectations as riboswitches with higher ARs typ- 
ically have lower gene expression in the absence of theo- 
phylline (22,28). Consistent with the observations of Figure 
2B, the riboswitches that display the highest ARs when a 
DNA template is used, those that we would predict to have a 
dominant kinetic component, show minimal theophylline- 
dependent increases in gene expression when an mRNA 
template is used in extract (Table 1, Supplementary Figure 
S7). However, the riboswitches with lower ARs (and less sta- 
ble OFF states), which we would predict to have a larger 
thermodynamic component, have comparatively higher lev- 
els of theophylline-dependent increases in gene expression 
when an mRNA template is used (Supplementary Figure 
S7), consistent with having a larger thermodynamic com- 
ponent. 

Finally, we used the opportunity of working in an in 
vitro system to address a long observed discrepancy between 
the reported of ~400 nM of the theophylline aptamer 
for its ligand (29,30) and the concentration of theophylline 
required to activate the riboswitch in intact E. coli cells 
(^300-500 fxM) (11,31). We performed dose-response re- 
lationships using the C-switch and a no-aptamer luciferase 
control both in intact E. coli cells and in extract. Surpris- 
ingly, the two dose-response curves obtained from extract 
and from cells are highly similar (Figure 5). The only sig- 
nificant difference occurs when the concentration of theo- 
phylline increases above 2 mM. In cells, luciferase activity 
begins to level off at 2 mM theophylline. However, in ex- 
tract, luciferase activity increases until at least 6 mM theo- 
phylline, a concentration that greatly reduces the extract's 
ability to produce functional luciferase, as seen in the dose- 
response relationship for the no-aptamer control (Figure 
5B). We and others (11,12) had previously interpreted the 
large difference between the of the aptamer for theo- 
phylline and the concentration of theophylline required to 
elicit changes in gene expression in cells resulted from a lack 
of permeability of the cells to theophylline (or an active ef- 
flux mechanism). In our in vitro experiments, theophylline 
permeability is not an issue, yet we still observe that mil- 
limolar concentrations of theophylline are required to ac- 
tivate the riboswitches. This result is consistent with a ki- 
netic trapping mechanism for riboswitch function, where 
the concentration of ligand at the time of initial folding de- 
termines the structure formed, not the amount of ligand 
present throughout the lifetime of the molecule. 

DISCUSSION 

We began this study to determine whether theophylline- 
dependent synthetic riboswitches function in vitro in the 
same manner that they function in intact cells. We discov- 
ered that the function of riboswitches encoded as DNA 
templates and expressed in E. coli S30 extract closely mim- 
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Figure 3. Decoupling transcription and translation in E. coli S30 extract. (A) Luciferase activity of the D riboswitch added to the S30 extract as a DNA 
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ics the function of these same riboswitches when expressed 
in intact E. coli cells. However, when these riboswitches 
were produced as mRNA templates by run-off transcrip- 
tion and then added to E. coli extract, they failed to show 
strong theophylline-dependent increases in gene expression. 
Our study of riboswitch function in extract demonstrates 
that the system serves as a good model for riboswitch 
function in vivo, can be used to gain a greater mechanis- 
tic understanding of synthetic riboswitches, and can likely 
be adapted to engineer novel riboswitches in a function- 
ally relevant environment. Our results also suggest that in 
vitro, these riboswitches operate by a kinetic mechanism, 
in which the 'decision' to translate an mRNA is made co- 
transcriptionally. 

These results were surprising to us because previous stud- 
ies from our lab (22) suggested that riboswitches from this 



family were under thermodynamic (equilibrium) control. In 
those experiments, we showed that a reporter gene under the 
control of F switch in intact cells produced theophylline- 
dependent changes in gene expression in the absence of ac- 
tive transcription, although these changes were relatively 
small. It is important to note that the functional properties 
of F switch are not affected by the reporter gene sequence 
(Figure 2). The data shown in Figures 2 and 3 show the 
existence of an important kinetic component to riboswitch 
function, suggesting that our earlier studies did not tell the 
entire story. 

Several of our previous observations of riboswitch func- 
tion make better sense when viewed in light of these 
new data. For example, we previously observed that nu- 
clease digestion experiments of theophylhne-dependent ri- 
boswitches showed only modest ligand dependent changes 
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Table 1. D family of riboswitches: names, variable sequence, free energy of the 5' UTR (kcal/mol), free energy of the sequester helix (SH, kcal/mol), and 
observed activation ratios either in vivo or as an RNA template in extract (Supplementary Figure S7) 



Construct 


Sequence 


-AG 5' UTR 


-AGSH 


In vivo AR 


In vzYroRNA AR 


615 


CACUUCUA 
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Figure 4. Plot of the observed activation ratios versus the m/o/tZ-predicted 
AG of the OFF state sequester helices for a set of 'D family' riboswitches. 
Each data point represents a single riboswitch sequence, including the par- 
ent D switch (open circle). The activation ratios were estimated using ex- 
pression data that can be found in Supplementary Table SI. 



in RNA structure (28). We had suggested that the modest 
changes resuhed from a small fraction of functional RNA 
sequences appearing against the background of a larger 
population of non-functional RNAs. We now argue that 
kinetic factors, which would not be detectable in experi- 
ments lacking active transcription (e.g. SHAPE; Tl nucle- 
ase mapping), contribute greatly to the function of these ri- 
boswitches and that this explains why only modest changes 
were previously observed. The current data suggest that the 
equilibrium component is only a small contributor to the 
overall function of these riboswitches. 

The present data also help explain why theophylline- 
dependent riboswitches expressed in cells typically respond 
to millimolar concentrations of ligand, while the aptamer 
binds theophylline at submicromolar concentrations in vitro 
(29,30). We previously attributed this difference to obser- 
vations regarding theophylline metabolism and permeabil- 
ity in E. coli (32) that suggested that the intracellular con- 
centration of theophylline relative to the medium may be 
~ 1000-fold lower (11). After correcting for this difference, 
it appeared that our earliest reported riboswitches acted 
under equilibrium conditions (11,22,31). However, several 
of our riboswitches display ARs that are inconsistent with 
a mechanism in which equilibrium phenomena dominate 
(26,28). In our in vitro experiments, the riboswitches show 
half maximal activation at nearly millimolar concentrations 
of theophylline. As ligand permeabihty is no longer an is- 
sue, these data further implicate a kinetic component to the 
mechanism. 



Our data suggest that this family of theophylline- 
dependent synthetic riboswitches operates via a mechanism 
that includes both a kinetic component (which requires ac- 
tive transcription) and an equilibrium component (which 
doesn't require active transcription). However, the overall 
behavior of the riboswitches is dominated by the kinetic 
component. Our results are consistent with a recent study 
by Forster et al that explored the importance of the kinet- 
ics of ligand binding in engineered tetracycline riboswitches 
(33). 

Most synthetic riboswitches have been created from ap- 
tamers that were selected for their ability to bind a ligand 
tightly. It has also been speculated that natural riboswitches 
may be derived from aptamers that were selected for their 
ability to bind ligands tightly in an RNA world, but were 
later repurposed to function in riboswitches. Consistent 
with this suggestion, the aptamers from the flavin mononu- 
cleotide (FMN) (34) and cyclic-di-GMP riboswitches bind 
their ligands with equilibrium dissociation constants that 
are orders of magnitude lower than the ligand concentra- 
tions that the riboswitches sense in cells. The cyclic-di-GMP 
aptamer, for example, binds its Hgand with a T^d of ~10 pM. 
It achieves such tight binding because its dissociation rate 
is extremely slow {t\/2 ~ 43 days) (35). As E. coli grown in 
rich medium divide approximately every 20 min and the av- 
erage mRNA half-life is ~6 min (36), this riboswitch can- 
not reach equilibrium, and thus acts as a kinetic trap like 
the theophylline riboswitches described here. 

These observations suggest that when there is a strong 
kinetic component to riboswitch activity, the equilibrium 
dissociation constant of the aptamer may not accurately re- 
flect riboswitch performance. In such cases, the Hgand as- 
sociation kinetics may be a better predictor of riboswitch 
function. This has imphcations for aptamer selection and ri- 
boswitch design. For example, an aptamer may display tight 
binding that is the result of a very slow ligand dissociation 
rate. However, such an aptamer may perform very poorly in 
the context of a kinetically controlled riboswitch. We sug- 
gest that aptamers that maximize ligand association kinetics 
are better suited for kinetically controlled riboswitches. 

The synthetic riboswitches reported thus far respond to 
only a small number of ligands. One of the biggest hurdles 
in the creation of synthetic riboswitches has been the dis- 
covery of new aptamers through in vitro selection that are 
appropriate for riboswitch function (2). Another hurdle is 
that aptamer selections are but one step in the process of 
discovering new riboswitches. Historically, aptamer selec- 
tions have been carried out in cell-free systems that allow 
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Figure 5. Dose response curves for luciferase reporters. (A) Dose-response curves for luciferase reporters expressed in intact E. coli cells. The gray squares 
correspond to the C riboswitch; the black squares correspond to a control lacking the aptamer. Data are normalized to the expression of either the C 
riboswitch in 10 mM theophylline or the control in 0 mM theophylline. Each data point represents the average of three independent measurements; error 
± SD. (B) Dose-response curves for luciferase reporters expressed from DNA templates expressed in E. coli S30 extract. The gray squares correspond to 
the C riboswitch; the black squares correspond to a control lacking the aptamer. Data are normalized to the expression of either the C riboswitch in 6 
mM theophylline or the control in 0 mM theophylline. The dose-response curve was terminated at 6 mM theophylline in extract because of the increasing 
toxicity of the ligand, as seen with the control dose response (black squares). Each data point represents the average of three independent measurements; 
error ± SD. 



one to assay ~10^^ sequences, while screens and selections 
for riboswitch function have been carried out in cells, which 
limits library sizes to ~10^ members. Our results and others 
(24,25) suggest that an in vitro system utilizing E. coli S30 
extract recapitulates riboswitch function in intact cells. This 
opens the door to directly selecting new riboswitches from 
large (^^10^^ members) libraries of RNA without perform- 
ing a separate aptamer selection. We expect that it will soon 
be possible to directly select riboswitches using in vitro sys- 
tems with the expectation that these riboswitches will read- 
ily function in intact cells. 
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